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Abstract. We have searched the Fourth Catalogue of 
Nearby Stars for halo stars and identified 15 subdwarfs 
and a high velocity white dwarf in the solar neighbour- 
hood. This search was motivated mainly by the recent de- 

\ terminations of MACHO masses of about 0.5A4q, which 
are typical for halo dwarfs. The local mass density of these 

I stars is 1.5-10 -4 A4q pc~ 3 , which is only 3% of the cur- 
rent estimate of the local mass density of the MACHO 
population. We compare the local density of subdwarfs 
with constraints set by HST observations of distant red 
dwarfs. Using models of the stellar halo with density laws 
that fall off like r~ a , a = 3.5 to 4, we find that the HST 
constraints can only be matched, if we assume that the 
stellar halo is flattened with an axial ratio of about 0.6. 
The non-detection of the analogues of MACHOs in the 
solar neighbourhood allows to set an upper limit to the 
luminosity of MACHOs of Mb > 21 magnitudes. 

Key words: Galaxy: halo - Galaxy: solar neighbourhood 
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1. Introduction 

The stellar halo of the Galaxy has been studied with re- 
newed interest since the MACHO (Alcock et al. 1996) 
and EROS (Ansari et al. 1996) collaborations have re- 
ported results of their campaigns of observing micro- 
lensing events towards the LMC. These indicate that MA- 
CHOS have masses typically of half a solar mass and con- 
tribute about one half to the mass budget of the halo of 
the Milky Way. Such masses are typical for red and white 
dwarfs, implying that MACHOs might possibly be mem- 
bers of the stellar halo instead of the dark halo of the 
Galaxy, as was previously thought. Bahcall and collabo- 
rators (Bahcall et al. 1994, Flynn, Gould & Bahcall 1996, 
Gould, Bahcall & Flynn 1997) have made by very deep 
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star counts based on HST data in situ measurements of 
the space density of such stars in regions, where the micro- 
lenses are expected to be physically located. Their conclu- 
sion is that halo dwarfs contribute only a few percent to 
the mass of the halo and are thus unlikely MACHO candi- 
dates. Graff & Freese (1996) infer from the same data that 
the local mass density of halo dwarfs is even less than one 
percent of the combined local densities of the dark and 
stellar halos. 

On the other hand, halo dwarfs are directly observed 
in the solar neighbourhood. Liebert (1995) discusses the 
results of the USNO parallax programme for red dwarfs 
and finds that the number of halo dwarfs in that sam- 
ple is consistent with the HST data. We have followed a 
complementary approach. Recently the Third Catalogue 
of Nearby Stars (CNS3) has been completed at the As- 
tronomisches Rechen-Institut (Jahreifi & Gliese 1997) and, 
after HIPPARCOS (ESA 1997) parallaxes and proper mo- 
tions have become available, is presently updated to the 
Fourth Catalogue of Nearby Stars (CNS4; JahreiB & Wie- 
len 1997). This provides the now most complete inventory 
of the solar neighbourhood up to a distance of 25 pc from 
the Sun and allows a new determination of the local den- 
sity of halo stars. Preliminary results based on the CNS3 
have been reported by Fuchs & Jahreifi (1997). 

2. Subdwarfs in the CNS4 

We have searched the CNS4 for halo dwarfs. These reveal 
themselves by their position on the subdwarf sequence in 
the colour-magnitude diagram (CMD), their low metallic- 
ity, and their high space velocities. Unfortunately, there 
is in this metallicity range a broad overlap of the com- 
paratively few halo stars and the much more abundant 
thick disk stars. We have therefore adopted a very con- 
servative search strategy in order to isolate a true halo 
population, even if this will lead to an undersampling of 
the halo population. We have examined first the stars on 
the subdwarf sequence, which can be clearly seen in the 
CMD shown in Fig. 1. The reliability of the photometry 
and the parallaxes of each star has been checked carefully, 
which led to the omission of quite a number of stars with 
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Fig. 1. Colour-Magnitude diagram 
of nearby stars with <tm v < 0.3 
mag. 



inaccurate parallaxes. In this way we identified LHS 2815, 
G1191, G1781, WO 9722 = LHS 64, GJ 1062, LHS 3409, 
and LHS 375 as likely halo stars. A number of stars in the 
CNS4 are classified as subdwarfs by their spectral type. 
This enables us to search for subdwarfs in that part of the 
CMD, V — I < 1.5, where the subdwarf sequence comes 
very close to the main sequence. We include GJ 1064A+B 
and Gl 53A+B into our list on the basis of this criterion. 
All stars, with the exception of LHS 2815, have decidedly 
large space velocities (cf. Table 1 and Fig. 2). In addition 
to these stars we found 5 stars, G1158, G1299, G1451A, 
WO 9371 = LHS 42 and the white dwarf Gl 699.1 (spectral 
type DA7) , which do not lie very prominently on the sub- 
dwarf sequence, but have space velocities, which clearly 
identify them as halo stars. 



All the stars in our list have very low metallicities. 
Leggett (1992) has determined the metallicities of a large 
number of nearby M dwarfs from their position in the 
(J — H) — (H — K) infrared-colour diagram. Accord- 
ing to that determination G1191, G1781, GJ1062, and 
G1299 have metallicities [Fe/H] < -1 dex. LHS 375 has 
been classified as a cold subdwarf by Ruiz and Anguita 
(1993). The metallicities of G153A, G1158, G1451, and 
GJ1064A, B have been taken from the compilation by 
Taylor (1994). Reid et al. (1995) have classified GJ 1062, 



G1781 0, WO 9722, and LHS 3409 as likely subdwarfs on 
the basis of their spectra. Similarly, Gizis (1997) has clas- 
sified WO 9371 as a subdwarf. Jones et al. (1996) con- 
firm the low metallicity of G1299. The metallicity of 
LHS 2815 has been determined by Axer, Fuhrmann & 
Gehren (1994). There are many apparently low metallic- 
ity stars in the catalogue, which must be thick disk stars. 
For instance Leggett (1992) has assigned to 20 stars in the 
CNS4 a metallicity of [Fe/H] < -1 dex, of which only 4 
have been finally included in our sample. 

In Table 1 we give the list of stars selected out of the 
CNS4 as described above, of which we are now reason- 
ably certain that they are genuine halo stars. The velocity 
components of the subdwarfs given in columns 6 to 8 in 
Table 1 are referred to the Sun. In order to reduce them 
to the LSR the solar motion (U, V, W) Q = (+9, +12, +7) 
km/s has to be added. Note that U points towards the 
galactic center. The velocity distribution of our sample is 
shown as a Toomre diagram in Fig. 2. It is consistent with 
a distribution centered on V = -220 km/s. No clustering 
near V = -40 km/s, which would indicate a contamina- 
tion by thick disk stars, is detected. The dispersion of 
y/U 2 + W 2 is 155 ± 30 km/s, which is within statistical 
errors consistent with other determinations (Layden 1995) 



1 G1781 is a variable star and shows H Q emission. However, 
it is a spectroscopic binary, which, in our view, accounts for 
this phenomenon. 
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Table 1. Subdwarfs in the CNS4. 
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HIPPARCOS parallax and proper motions 



of the kinematical parameters of halo stars. The masses cally for a metallicity of [Fc/H] = -1.35 dex by Alexander 

et al. (1997). 

The changes due to the introduction of HIPPARCOS 
parallaxes can be seen by comparing Table 1 and the cor- 
responding Table 1 of Fuchs & Jahreifi (1997). 

3. Results and discussion 

3.1. Local density 

The local mass density of halo stars can be estimated in 
various ways. In Fig. 3 we show density estimates, which 
have been calculated from the cumulative distribution of 
the stars in our sample, i.e. by adding up the masses of 
stars within given distances and dividing by the corre- 
sponding spherical volumes. The stars within 10 pc give 
rather high density estimates of the order of 1CP 3 A4q 
pc~ 3 . Wielen & Jahreifi (Wielen 1976) found similar val- 
ues, when evaluating the second edition of the Gliese cat- 
alogue. A careful analysis of the spatial distibution of the 
stars shows, however, that the high density peaks are of 
purely statistical nature and certainly not realistic. We 
find a plateau in the density run, when we consider the 
stars within the 20 pc sphere (cf. Fig. 3). Beyond that 
the density drops off, probably because the catalogue be- 
comes incomplete. The most likely value of the local mass 
density lies according to our determination in the range 
1.5 to 1-10~ 4 Mq pc~ 3 . We note that just outside the 
nominal boundary of the CNS4 there are two further sub- 
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Fig. 2. Toomre diagram of the nearby halo dwarfs. The dashed 
lines indicate lines of constant kinetic energy of the halo stars. 



given in the last column of Table 1 have been assigned 
using the mass-to-luminosity relation calculated thcoreti- 
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dwarfs, LHS 173 and GL 788.2, at distances 25.5 pc and 
25.6 pc, respectively. If they are included, the mass density 
is 1-10 - Mq pc -3 , which, in our view, is a firm lower limit 
to the local mass density of halo stars.. The local density 
of the dark halo (Bahcall, Schmidt & Soneira 1983, Gates 
et al. 1996) is about 9-1CT 3 M Q pc~ 3 , so that the local 
density of halo stars is about 1.7% of the halo density or 
3% of the local MACHO density as determined by the 
MACHO collaboration. We have broken down our sample 
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Fig. 3. Mass density of nearby halo dwarfs as function of the 
distance from the Sun. Within 10 pc for some of the stars 
only the positions are indicated. The errorbars indicate the 
statistical uncertainties. 



of subdwarfs with respect to absolute magnitudes. The 
resulting 'luminosity function' is compared with the sub- 
dwarf luminosity function obtained by Dahn ct al. (1995) 
for stars with My > 10 mag in Table 2. Both agree well 
within statistical uncertainties. 



Table 2. Luminosity function of red subdwarfs. 
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3.2. Matching the constraints by HST data. 

It is interesting to confront the density of halo stars de- 
rived here with density estimates based on the HST data. 
Graff & Freese (1996) find in a detailed analysis of earlier 
HST data (Bahcall et al. 1994) a value about five times 
less than our estimate, which is partly due to their restric- 
tion to red dwarfs redder than V — I ~ 2 mag. Whereas 
these authors tried to infer the local density of halo stars 
from the star counts, we reverse the approach and make 
predictions of the expected number of stars in the star 
count field. Recently Flynn et al. (1996) have reported an 
analysis of the Hubble Deep Field (HDF). The HDF has a 
size of 4.4 square arcminutes and is located at I — 126°, b 
= 55°. Flynn et al. (1996) were able to identify stars down 
to a limiting magnitude of 26.3 mag in the /-band. 

In the following we assume for the population of halo 
stars a density law of the form 



r 2 _|_ r 2 



(1) 



with an arbitrarily chosen core radius of r c =1 kpc. 
Such density laws are typical for other tracers of the halo 
population, such as RR Lyrae stars (Wetterer & McGraw 
1983), horizontal branch stars (Preston et al. 1983), in the 
Milky Way or external galaxies (Pritchet & van den Bergh 
1994). The exponent a lies typically in the range 3.5 to 4. 
Note that the dark halo is usually described by a density 
law of the form as in equation (4) with a = 2. If the sub- 
dwarfs are observed locally at a density ^h Q , one predicts 
by integrating along the line of sight 



■n 



Jd min V 



(2) 



r c + r G) + — 2sr Q cos I cos b 



ds 



stars in the star count field. f2 is the angular area of the 
field and the minimum and maximum distances are deter- 
mined by the limiting magnitudes. 

We concentrate first on the colour range V — I = 4.8 
to 3.5 mag. We have 4 stars in that range within 16.5 pc, 
where we believe our sample to be reasonably complete. If 
we project these into the cone towards the HDF, we obtain 
the star numbers summarized in Table 3. The limiting 
distances have been calculated for each star individually, 



d 



mm, max 



10 0.2(J 



mm, max 



+(y-/)-M v +5) 



(3) 



where, in order to avoid confusion with disk stars in the 
HDF, we adopt a lower limit of I m i n = 24.6 mag. As can 
be seen from Table 3 we predict 7 to 11 stars in the HDF, 
whereas Flynn et al. (1996) have actually detected no star. 
Despite the low number of stars involved, this discrepancy 
seems to be statistically significant. Several explanations 
might account for this discrepancy. First, the stellar halo 
might be much more irregular and lumpy than previously 
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Table 3. Predicted Number of Stars in the HDF and the GS. 



V -I > 1.8 
a c/a riHDF 


V-K1.8 

riHDF nGS 


3.5 1 11 
4 1 7 

3.5 0.6 4 
4 0.6 2 


17 258 
7 140 
5 67 
2 31 



thought. Second, the stellar halo is almost certainly flat- 
tened (Wetterer & McGraw 1983, Preston, Shcctman & 
Beers 1983, Pritchet & van den Bergh 1994) with an axial 
ratio around c/a w 0.6. If we take such a flattening into 
account in the halo model, 



+ R 2 r. 



-l+R 2 + z 2 /{c/a) 2 



(4) 



where R, z denote cylindric coordinates, we predict 2 to 4 
stars in the HDF, i. e. star numbers which are statistically 
consistent with no star seen by Flynn et al. (1996). 

The star counts in the HDF in the colour range V — I < 
1.8 mag can be interpreted in a similar way. Disk stars in 
this colour range are so bright that they would have to 
lie several kpc above the midplane to appear fainter than 
I = 22 mag. Thus in order to avoid confusion with disk 
stars in the HDF, we consider only stars fainter than I = 
22 mag. In our sample there are 4 stars within 20 pc in 
this colour range. Their V — I colours cluster around 0.9 
and the white dwarf has V — I = 0.5. Thus we define a 
colour range 0.5 < V — I < 1.8, in which we compare the 
predictions with actually observed numbers of stars. If the 
stars of our sample are projected into the cone towards the 
HDF using spherical halo models, we predict 7 to 17 stars 
depending on the model, whereas Flynn at al. (1996) have 
observed 4 stars. If we assume again a flattening of c/a = 
0.6, the number of expected stars is consistent with the 
observed number. 

Gould et al. (1997) have very recently investigated an- 
other field observed by HST. This is the so-called Groth 
strip (GS), which has an angular size of 114 square arcmin- 
utes and is located at I = 96° and b = 60°. Its limiting 
magnitude is I = 23.9 mag. Gould et al. (1997) have very 
kindly made their data available to us, so that we were 
able to make star counts in this field. Due to the limiting 
magnitude of I = 23.9 mag even the faintest stars in the 
GS with colours redder than V — I = 1.8 could be disk 
stars. Thus we focus on the colour range V — I < 1.8 and 
adopt the same halo zone as in the HDF. The GS has 66 
stars in this zone. There is a slight chance that this zone 
might be contaminated by faint white dwarfs in the galac- 
tic disk. It can be shown, however, that, if the local density 



of disk white dwarfs as determined by Jahreifi (1987, 1997) 
and conventional estimates of the vertical scale height of 
the galactic disk are adopted, one would expect about one 
white dwarf in this zone of the GS. If we project the local 
subdwarfs into the cone towards the GS, we obtain the 
numbers of predicted stars summarized in the last column 
of Table 3. Again a flattened halo model is a better fit to 
the data. 

3.3. Constraints on the absolute magnitude of MACHOs. 

Like previous authors we have not detected the analogues 
to MACHOs in the solar neighbourhood. This raises the 
question why. All subdwarfs in our sample appear in the 
LHS high-proper motion star catalogue. Being halo ob- 
jects, MACHOs in the solar neighbourhood must have 
similar proper motions. The fact that they have not been 
found in high-proper motion surveys allows to set an up- 
per limit to their luminosity (Fuchs & Jahreifi 1997). The 
LHS catalogue (Luyten & La Bonte 1973, Dawson 1986) 
is claimed to be about 90 % complete for stars down to 
apparent magnitude B — 21 mag and with proper motions 
2. 5/a > fi > 0. 5/a. For stars brighter than B = 10 mag 
it is claimed that all stars with proper motions fi > 0. 3/a 
are contained in Luyten's NLTT catalogue. We assume 
that the MACHOs are homogeneously distributed around 
the Sun and have a gaussian velocity distribution, 



dN 



VMr. 



(2^) 3 /V : 



M 



1 



(U 2 +(V-V) 2 +W 2 )d 3 vd 3 r ,(5) 



centered on V = -220 km/s and with a velocity disper- 
sion typical for halo objects, ctm = V / \/2- Integrating now 
over all radial velocities and the proper motions accord- 
ing to the specifications of the LHS we can determine the 
radius r 2 of the sphere out of which one would expect say 
2 MACHOs in the LHS, 



4.74 
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2^0 / J m „4 
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°m Jo Jo 



■To ( -^-4.74/xrVl - sin 2 I cos 2 b ) 

■exp - -V (V 2 (l - sin 2 I cos 2 b) + (4.74^r) 2 ) , (6) 



where I denotes the Bessel function. According to Poisson 
statistics two expected MACHOs would be still consistent 
with no MACHO actually seen in the LHS survey. In Ta- 
ble 4, we give results of numerical integrations of equation 
(6) assuming a local density of MACHOs of vm & = 0.5 • 
0.01 M Q pc~ 3 / 0.5 M Q as determined by the MACHO 
collaboration. A lower limit of the absolute magnitude of 
the MACHOs is then given by Mb = 21 - 51ogr 2 + 5, 
because MACHOs more distant than r 2 will not show up 
in the LHS, since they are fainter than the apparent mag- 
nitude limit of the LHS, B = 21 mag. From Table 4 we 
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Table 4. Estimated absolute magnitudes of MACHOs. 
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20.1 


18 
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conclude that MACHOs are fainter than Mb = 21.2 mag. 
The southern hemisphere (5 < —30°) and the galactic belt 
are not sampled by the LHS survey. Thus one could argue 
that only the region b > 30° is properly sampled. This 
would reduce the numbers in the second column of Table 
4 by a factor of 0.25 and shift the lower limit of the ab- 
solute magnitude of the MACHOs to M B = 20.6 mag. If 
the distribution of MACHOs is irregular and lumpy, the 
micro-lcnsing results could mimic a too high mean den- 
sity of MACHOs. If we assume a local MACHO density 
of only 10% of the value determined by the MACHO col- 
laboration, the lower limit of the absolute magnitude of 
MACHOs would be M B = 20.1 mag. 

If the MACHOs were identified with faint red dwarfs, 
their masses would have to be less than 0.1 A4® in con- 
tradiction to the estimate by the MACHO collaboration. 
As an alternative very faint white dwarfs have been dis- 
cussed in the literature as MACHO candidates (Chariot & 
Silk 1995). Graff et al. 1997 discuss various aspects of this 
scenario and point in particular to the problem that such 
faint white dwarfs must have been born due to the long 
cooling times very early in the evolution of the universe. 

What would be the requirements to detect such faint 
objects in future proper motion surveys? Assuming, for 
instance, a limiting magnitude of B = 22.5 mag as for 
UKST plates, an object with an absolute magnitude of 
Mb = 21 mag could be detected up to a distance of 20 pc 
from the Sun. Its typical proper motion would be 3. 5/a. 
Assuming again a volume density of 10~ 2 pc~ 3 , one would 
expect 333 objects distributed evenly over the sky. 

Acknowledgements. We are indebted to C. Flynn, K. Freese, 
N.W. Evans, and R. Wielen for many valuable hints and com- 
ments. We are grateful to J. Bahcall, C. Flynn, and A. Gould 
for making available to us the their data on the Groth strip 
prior to publication. This research has made use of the SIM- 
BAD database, operated at CDS, Strasbourg, France. 

References 



Axer M., Fuhrmann K., Gehren T., 1994, A&A 291, 895 
Bahcall J.N., Flynn G, Gould A., Kirhakos S., 1994, ApJ 435, 
L51 

Bahcall J.N., Schmidt M., Soneira R.M., 1983, ApJ 265, 730 
Chariot, S., Silk J., 1995, ApJ 445, 124 

Dahn GC, Liebert J, Harris H.G, Guetter H.H., 1995, in: 'The 
Bottom of the Main Sequence and Beyond', ed. GG. Tinney 
, Springer, Berlin, 239 

Dawson P.G, 1986, ApJ 311, 984 

ESA, 1997, 'The HIPPARCOS Catalogue', ESA SP-1200, ES- 
TEC, Nordwijk 

Flynn G, Gould A., Bahcall J.N., 1996, ApJ 466, L55 

Fuchs B., Jahreifi H., 1997, in: 'Aspects of Dark Mat- 
ter in Astro- and Particle Physics', eds. H.V. Klapdor- 
Kleingrothaus, Y. Ramachers, World Scientific, Singapore, 
in press 

Gizis J.E., 1997, AJ 113, 806 

Gould A., Flynn G, Bahcall J.N., 1997, ApJ, submitted 
Graff D., Freese K., 1996, ApJ 456, L49 

Graff D., Laughlin G, Freese K., 1997, ApJ, submitted/ as- 

troph 9704125 
Jahreifi H., 1987, Mem.S.A.It. 58, 53 

Jahreifi H., Wielen R., 1997, in: 'Presentation of the HIPPAR- 
COS and TYCHO catalogues', eds. M.A.G Perryman, P.L. 
Bernacca, ESA-SP 402, ESTEC, Nordwijk 

Jahreifi H., Gliese W., 1997, in preparation; preliminary ver- 
sion available on CD-ROM 'Selected Astronomical Catalogs 
Vol. 1', eds. L.E. Brotzman and S.E. Gessner, ADC, NASA, 
Greenbelt 1991 

Jones, H.R.A, Longmore A.J, Allard F., Hauschildt P.H., 1996, 

MNRAS 280, 77 
Layden A.G, 1995, AJ 110, 2288 
Leggett S.K., 1992, ApJS 82, 351 

Liebert J., 1995, in: 'Dark Matter', AIP Conf. Proc. 336, 71 
Luyten W.J., La Bonte A.E., 1973, The South Galactic Pole, 

Univ. of Minnesota, Minneapolis 
Preston G.W., Shectman S.A., Beers T.G, 1983, ApJ 375, 121 
Pritchet C.J., van den Bergh S., 1994, AJ 107, 1730 
Reid I.N., Hawley S.L., Gizis J.E., 1995, AJ 110, 1838 
Ruiz M.T., Anguita G, 1993, AJ 105, 614 
Taylor B.J., 1994, PASP 106, 704 
Wetterer C.J., McGraw J.T., 1983, A J 112, 1046 
Wielen R., 1976, Mitt. Astron. Ges. 38, 254 



Alcock G et al., 1996, ApJ 461, 84 
Alexander D.R. et al., 1997, A&A 317, 90 
Ansari R. et al., 1996, A&A 314, 94 



This article was processed by the author using Springer- Verlag 
DT E X A&A style file L-AA version 3. 



